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We consider resonant dynamics in a dilute atomic gas falling under gravity through a periodically pulsed 
standing-wave laser field. Our numerical calculations are based on a Monte Carlo method for an incoherent 
mixture of noninteracting plane waves, and show that quantum resonances are highly sensitive to the relative 
acceleration between the atomic gas and the pulsed optical standing wave. For particular values of the atomic 
acceleration, we observe fractional resonances. We investigate the eff'ect of the initial atomic momentum width 
on the fractional resonances, and quantify the sensitivity of fractional resonances to thermal eff'ects. 
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I. INTRODUCTION 

The atom-optical kicked rotor can be realized by subjecting 
a dilute atomic gas to a pulsed optical standing wave. This 
system has become well established as a convenient way to 
study quantum chaotic dynamics |[Il[2||3l|4l[5l[6|l7,8,9,10, 

[iTifnifTsiiiiisiiisiiniiiaiisEa^ 

[28, 29, 3Ql [311 [321 [331 El- In particular, quantum resonances 
and antiresonances ES [36l |37l|38| [39l |W, 41 42, 43,^113 
are a dramatic signature of quantum chaos [46^47^481, and 
these phenomena have been observed and studied in detail in 
cold atom kicked-rotor experiments d 121 [3 |4]|. One of the 
features of such experiments is that the optical standing wave 
is often oriented vertically, allowing the atom cloud to fall un- 
der gravity during the optical pulses. This realizes the quan- 
tum kicked accelerator |4|, and has led to the observation of 
quantum accelerator modes |[4|[5l[6l|7||8l|9|[T0l[II|[ia which 
are closely related to quantum resonances. However, quantum 
resonances themselves have not been specifically investigated 
for the quantum delta-kicked accelerator. 

In this paper, we consider the atom-optical quantum kicked 
accelerator where the linear time-independent potential, typi- 
cally provided by the gravitational acceleration, can be freely 
chosen. This is possible in experiments with either verti- 
cal or horizontal orientation of the optical standing wave, as 
an appropriate variation in the phase of the optical standing 
wave can induce an eflTective gravitational acceleration. We 
find that the relative acceleration between the atom cloud and 
the pulsed optical standing wave can lead to fractional reso- 
nances. We characterize these fractional resonances and con- 
sider in detail the eflTect of the initial atomic momentum width 
on the fractional resonant dynamics. This generalizes our pre- 
vious discussion on the eflTect of temperature on quantum reso- 
nances and antiresonances in the atom-optical quantum kicked 
rotor ll49i] . 

This paper is organized as follows. In Sec. [ll| we derive 
the quantum ^-kicked accelerator Hamiltonian for the system 
of a two-level atom falling through an optical standing wave. 
We then consider the time evolution of the system based on a 
Floquet operator approach. Our treatment closely follows the 



derivation presented in Sec. II of Saunders et al. |49|. In Sec. 
Ill we consider fraction resonances in the zero-temperature 
limit and for the case of a broad, thermal atomic momentum 
distribution. As in our previous work 1 49 J , we use the evo- 
lution of the momentum moments to characterize the atomic 



dynamics. In Sec. [IV| we consider the dynamics of a plane 
wave state and investigate the eflTect of the initial momen- 
tum on the system behavior. From that discussion we develop 
an understanding of the temperature dependence of fractional 
resonances. We conclude in Sec.[Vl 



II. ATOMIC SYSTEM 
A. Model Hamiltonian 

1. Freely-falling two-level atom in a laser field 

We consider a cloud of trapped and laser-cooled thermal al- 
kali atoms. The cloud is released from all external fields, and 
then falls through a vertically aligned pulsed sinusoidal poten- 
tial ||4|[5l[6l[7l[8ir9, 10, 11 1 formed by two counter-propagating 
laser beams |50|. We neglect interatomic collisions, as is rea- 
sonable for a dilute thermal gas. For the case of an initially 
Bose-condensed sample where atom- atom interactions can be 
significant, it is possible to tune the scattering length to zero 
[ 51 1 near a Feshbach resonance 1521 [531(541 [55l . 

We describe the system using a one-dimensional model 
along the vertical axis. As in |49|, the Hamiltonian H de- 
scribes a single two-level atom of mass M, with internal 
ground state |^) and excited state \e), but now the centre-of- 
mass motion is influenced by a linear gravitational potential. 
The internal atomic levels are coupled by two equal-frequency 
laser fields with phases (pi and 02, respectively, and laser 
wave- vector magnitude in the vertical direction ki. We con- 
sider the lasers to be of equal intensity, and far detuned so that 
spontaneous emission can be neglected. As in B91 , we work 
in a rotating frame defined by Ui = Qxp(i[ajL\e}{e\-oJo(\e}{e\-\- 
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\g){g\)l2]t), and so 



3. Transformation to a spatially-periodic Hamiltonian 



2M 



Uaj, cosihz + (Pr/2) [e'^^'\e}{g\ + H.c] . 



(1) 



Here z is the vertical centre-of-mass atomic position operator 
with p its conjugate momentum, t is the time, g is the local 
gravitational acceleration, 0^ = 0^ - 02, (ps = (0i + 02)/2, 
and the detuning A = ojq-col, with fiojo the energy difference 
between the internal atomic levels and ojl the laser frequency. 



2. Transformation to an accelerating frame 

In contrast to f49\, we consider the laser phases 0i and 02 
to be time-dependent. In particular, we choose 0^ = so 
that the laser field forms a "walking wave" with acceleration 
-^0 ISOl . This allows for the possibility of exactly negating 
the effect of the gravitational force by taking = g, and for 
the versatility of investigating different effective values of the 
gravitational acceleration by choosing at our convenience. 
The relative phase 0^ can be tuned with an accuracy on the 
order of one part in a million fSOl, allowing for precise control 
over the effective gravitational acceleration |[9l. 

The Hamiltonian ([T]) can be transformed to the frame ac- 
celerating with the walking wave by implementing the unitary 
transformation U2 = exp(i[Ma(pzt - a(ppt^/2 -h ^(t)]/fi), where 
for convenience ^(t) = Ma(pt^(a(p - 2g)/l2. This yields 



H' =nA\e}{e\ + ^ + 



^U^RCOsihz) [e'^^\e){g\^n.c]. 



(2) 



where we have used U2ZUI 



a(pf/2 and U2PUI = p 



Ma^pt, and we have defined a = g - a^. It is evident from Eq. 
([2]) that, in the frame defined by U2, a assumes the role of an 
effective gravitational acceleration. 

In the far-detuned limit (Qr/A <$c 1), the excited state \e} 
can be adiabatically eliminated I421ISS1- By also invoking the 
unitary transformation U3 = exp(-/Q^|g)(g|r/8A), the Hamil- 
tonian simplifies to 



p2 f^Q2 



(3) 



where K = 2kL, and the atoms are assumed always to be in 
the internal ground state. In an exactly analogous fashion to 
f49 |, we consider the laser standing- wave amplitude to be pe- 
riodically pulsed, with period T and pulse duration tp. For 
sufficiently short pulses |[T3l[T71|49l, It is then possible to ap- 
proximate the Hamiltonian by 



^c5ka = + Maz - Hd cos(Kz) ^ S(t - nT\ (4) 



where we have defined an effective potential depth 0^ = 
Q^tp/SA. A system governed by Hamiltonian ^ is known 
as a ^-kicked accelerator t4J. 



It is convenient to transform Hamiltonian (|4]) to a frame in 
which it takes a spatially periodic form so that we can apply 
Bloch theory f49l[57||5t|[59l[60l. This is accomplished by a 
unitary transformation defined by the operator L58. .59. 60j 



1/4(1) = Qxp(iMazt/fi). 
The transformed Hamiltonian is then 
(p - Matf 



(5) 



2M 



fKpd cos(Kz) ^ S(t - nT). (6) 



We expand the position and momentum operators into discrete 
and continuous components, i.e.. 



z = K-\2nl + 0), p = nK(k 



(7) 



where the eigenvalues of / and k are integers, the eigenval- 
ues ofOe [-71, 7i) and the eigenvalues of the quasimomentum 
operator y§ g [-1/2, 1/2). Substituting the reformulated oper- 
ators into Hamiltonian ([6]) yields 



[fiKik Matf 
2M 



- ficpd cos(e) ^ S(t - nT). (8) 



The quasimomentum operator y§ commutes with Hamilto- 
nian ([8]) |49 , 58 1 . Therefore, the quasimomentum is conserved 
in the accelerating frame, and the laser field induces coupling 
only between momentum eigenstates differing in momenta by 
integer multiples of HK, i.e., integer multiples of two photon 
recoils 141. 



B. Time evolution 

1. Transformed Floquet operator in the accelerating frame 

For a time-periodic system, the unitary time-evolution op- 
erator for a single temporal period is known as the Floquet op- 
erator. Hamiltonian ^ has temporal period T and we choose 
to define the Floquet operator from just before one kick to just 
before the next, i.e.. 



F^expl-- 



r exp(/0jcos(i^z)). (9) 



In the accelerating frame defined by Eq. ([5]), the Hamilto- 
nian ([6]) is not periodic in time, so strictly speaking there is no 
Floquet operator. However, it is useful to consider the trans- 
formed Floquet operator which is the time-evolution operator 
describing kick-to-kick dynamics of the system in the acceler- 
ating frame. The Floquet operator ([9]), transformed according 
to the unitary operator ([5]), becomes 



F{nT, [n - \]T) = U4(nT)FUl([n - l]T). 



(10) 
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The time variables on the left-hand side are required to in- 
dicate explicitly that the transformed Floquet operator corre- 
sponds to the system evolution from just before the (n - l)th 
kick to just before the ^th kick [ 58 1. 



To evaluate the transformed Floquet operator ( ^Q) ), it is con- 
venient to separate the terms that depend on z from those that 
depend on p. The Floquet operator ^ can be rewritten as 



F = exp 



6^ 



exp 



2M 2 



iMazT 
X exp I — -h 14 



I cos(Kz) , 



(11) 



as shown in Appendix [A| Substituting Eq. ( pTj ) into Eq. (10), 



then using U4(nT)pUl(nT) = p - ManT and simplifying, we 
find that 



F(nT, [n-l]T)= exp ( [^n^ -3n-\-l]T 



X exp 



n 



2M 



f(2«-i)r^ 



X exp (i(pd cos(Kz)) . 



(12) 

the global 
we absorb 



The first exponential term in Eq. (12), i.e., 
phase, ^ is not generally of interest. Therefore 
it into the transformed Floquet operator by defining F„ = 
Qxp(iMa^[3n^ -?>n+ V[T^ l6n)F{nT, {n - l]r), such that 



Fyi = exp 



pa 



2M 



{In - l)r^ j exp {i^d cos(Kz)) . 

(13) 

Finally, we substitute the discrete and continuous compo- 
nents of z and p from Eq. ^ into Eq. ([13]). A further simplifi- 
cation is then possible due to quasimomentum conservation 
(see Sec. [irA3l ) |58l |59l |60]|. When the transformed Flo- 
quet operator acts on a general quasi-momentum eigenstate 
\^(J3)} = Tjk <^k\k + yS), where the dimensionless momentum 
eigenkets \k -\- P) satisfy p\k + P) = fiK{k + P)\k + P) and 
{k' -\-p'\k-\- P) = Skk'Sip - p') 1 49 1, the operator p can be 
replaced by its eigenvalue p. Hence, the transformed Floquet 
operator, restricted to acting on a particular quasimomentum 
subspace I58l , is 



Fn(P) =exp|-- 



2M 

X exp (i(f)d cos(Kz)) . 



ik + pfT - Tinaik + P){2n - 1) 



(14) 



In Eq. ( 14 ), we have defined Q = KaT^ lln. This parameter 
is a dimensionless rescaling of the eff'ective gravitation accel- 
eration a = g - (20, in the frame comoving with the laser walk- 
ing wave (see Sec. II A 2 ). In the context of quantum accelera- 
tor modes, D has been referred to as the unperturbed winding 



^ The global phase is corrected slightly from that given previously 14]. 



number fTOl. This is in analogy with the sine-circle map fSTl, 
which can be considered a particular dissipative limiting case 
of the dynamics of the quantum ^-kicked accelerator. 



2. Quantum resonances and antiresonances 



The time evolution described by the transformed Floquet 



operator ( 14 ) varies significantly depending on the system pa- 
rameters. In this paper, we investigate the eff'ect of the ef- 
fective gravity parameter Q on quantum resonant dynamics in 
the strong-driving regime (K^Tfi(pd/M > 1.5). In the ^-kicked 
particle system (i.e., the ^-kicked accelerator with D = 0), 
quantum resonances and antiresonances occur for pulse peri- 
odicities T equal to integer multiples of the half-Talbot time 
Tt/2 1621, i.e.. 



^ 2 



^2nM 



(15) 



where ^ is a positive integer (see B9l and references therein). 
In this paper we consider in detail pulse periodicities given by 
Eq. ( p?] ) only. 



The transformed Floquet operator ( 14 ), for the case where 
the pulse periodicity T is given by Eq. (15]), simplifies to 



r:[k+2kj3£-m(2n-l)] ^icp^ cos 6 



where, because the eigenvalues of k are integer, we have 
used that exp(-/M2) = exp(-/M) EH [59l [60l. In the 
case where T = £Tt/2, the eff'ective gravity parameter is 
Q = 2nfaM'^in^K^. 

The transformed Floquet operator ([16]) can be written in a 
more concise way by defining 



Kjn =T:[{\+2p)£-Q.{2n-\)\, 



I.e., 



FniP) 



J[(l+l3)n£-Ky,]l3-iKy„kJ(Pd cos 6 



(17) 



(18) 



The first term in the transformed Floquet operator (18) pro- 
vides a quasimomentum dependent phase and the third term 
describes the momentum kick due to a laser field pulse. Con- 
sidering the second term in Eq. ( [TS] ) it is evident that, for an 
integer value of Q, Qxp(-iKyrik) Qxp(-iKyok) and the n- 
dependence of the transformed Floquet operator ( [TS] ) drops 
out. Furthermore, when P = 0, Kjq 7i{£ -\- Q), and we find 
that integer changes in the value of Q can be equivalently de- 
scribed by modifying the number of half-Talbot times making 
up a kick period T, i.e., we can define an eff'ective kick pe- 
riod T = tTT/2 from the integer t = f -\- ^. Recall that, 
for the ^-kicked rotor where D = 0, even (odd) values of 
£ are associated with resonant (antiresonant) behaviour 1491. 
Therefore, we expect that for a sufficiently narrow atomic mo- 
mentum distribution centred around zero, quantum resonance 
(antiresonance) will be observed for even (odd) Indeed, a 
change in Q by an odd integer value will change the system 
dynamics from resonant to antiresonant (or vice versa). 
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3. Rational values ofQ 



formally, 



The main focus of this paper is on rational values (particu- 
larly non-integer values) of Q, i.e., fl = r/s, where r and s are 
integers with no common factors and 5- > 0. To illustrate why 
these particular values are of interest it is useful to make one 
further transformation of the Hamiltonian ([6]) into the frame 
where the atom cloud is on average stationary, i.e., we apply 
the unitary operator ^5 = exp(-/[a^^^/2 - Ma^P /6]/fi) to the 
Hamiltonian ([6]). The transformed Hamiltonian is 

Hsm = ^-^^d cos(Kz - n^n^) V S(t - nT). (19) 

n-O 

In this frame, we see that the effect of Q., when the ^th kick 
is applied, is to shift the phase of the standing-wave poten- 
tial by an amount n^nfl. For the case where = r/s, this 
introduces an additional temporal periodicity into the system 
because the phase of the optical standing-wave pulses repeats 
identically with a period denoted by ht kicks. The additional 
periodicity can equivalently be deduced from the Floquet op- 
erator ([18]). In Appendix[B]we show that ht = s(l-\-sr mod 2), 
i.e., Ht = s for even sr, and Ht = 2s for odd sr. The Hamilto- 
nian ^T9\ can then be rewritten as 



/v2 rij-l 

^stat = ^-^(Pd X cosiKz-n'^-nr, 
2M ^ 

n'=0 



/s)^S(t-[nTn-\-n']T), 
(20) 

which applies for all possible values of T. 

For the case where Q is an even integer, ht = I and the 
standing- wave pulses are always in phase. It follows that the 
atomic dynamics are unchanged from the D = case. When 



Q is an odd integer, ht = 2 and from Eq. ( 20 ) we find that con- 
secutive optical pulses are exactly tt radians out of phase. We 
primarily consider T = £Tt/2, in which case the dynamics for 
even Q will always be either resonant or antiresonant depend- 
ing on whether £ is even or odd, respectively. Changes in by 
an odd integer cause the dynamics to swap between resonant 
and antiresonant, as described in Sec. |IIB 2"l and addressed fur- 
ther in Sec. |III A"3l Fractional values of Q lead to fractional 
resonances and antiresonances. Fractional resonances will be 
the main focus of this paper (see Sec.[Tll|). 



4. Momentum eigenstate evolution 

It is convenient to first determine the time-evolution of the 
system when the centre-of-mass state of the atoms is initially 
prepared in the momentum eigenstate 



\^^{t = Q)) = \k+p). 



(21) 



Once the evolution of state ( 21 ) is known, the evolution of any 
general state can be determined, as illustrated previously for 
theQ = Ocase 1491 . 

At time t = nT, the evolution of state ( [2T] ) can be deter- 
mined by consecutively applying transformed Floquet oper- 
ators defined by Eq. ([T8|, for each of the n periods. More 



where 



mt = nT)} = rnmk+/i}, 



(22) 



(23) 



Due to quasimomentum conservation (see Sec. II A 3 (, Eq. 
^22\ can be expanded using momentum eigenstates witli tlie 
initial quasimomentum yS, i.e., 



mt = nT)}= £ Ckj(Ii,nT)\j + p). 



(24) 



where tlie probability amplitudes Ckjifi, nT) are given by 

Ckjdi, nT)6(J3 - p') = (j + p'\fnmk + P)- (25) 



Evaluation of the matrix elements ( [25] ) closely follows our ear- 
lier derivation presented for the D = case |49|, and is given 
in detail elsewhere |[63]| . In summary. 



Ckj(fi,nT) = Jj-k(oj)e' 



.i(j-k)x^-2inrkJn^7:(k+m^-in7Tfi^i 



where^ 



for 



n-l 



and 



y = i^-i(2nr-n^nQ) ^ ^i(2jr-fni 
j=0 



Y = -nil + 2fi)L 



(26) 



(27) 



(28) 



(29) 



We will often specify the parameter dependence of v by writ- 
ing y = yn,{(^,P). 



The sum in Eq. ( [28] ) can be evaluated analytically for par- 
ticular choices of the dimensionless eff'ective gravitational ac- 
celeration Q. For integer values of D, i.e., Q = ri 1631 , 



v«/(ri,^) = ie 



-/(«+l)(T-;rri/2) 



sin(^[T - 7rri/2]) 
sin(T-7rri/2) 



(30) 



In Appendix [C] we consider half-integer values of Q, i.e., Q = 
r2/2. We find that 

ynAr2/2,P) = 2/.--(^t-i^--/4)cos(T-;rr,/4)5!^^ (31) 

sm(2T) 



Note that Eq. (31) applies for even values of n only. 



^ The parameter y in this paper is related to a previously defined parameter // 
[49] according to y = ijj.* . 
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In theyS = subspace, where T 7t£/2, v can be evaluated 
analytically for Q = 1/^, at kick values n which are specific 
multiples of s | 63 1. For ^(1 - £) even, and n a multiple of s, 



yn,Kl/s,0) = ie 



_ •-i7:[ni-n^Q.-{l-se^)IA] 



(32) 



For ^(1 - £) even or odd, and n a multiple of Is, 



Vn,i{Ms,0) = ie 



2^fs 



(33) 



From Eq. ( [33] ) two main classes of behaviour can be identi- 
fied. For odd values of ^(1 - i), Eq. (33 ) collapses to zero, 
i.e., the initial state is reconstructed every 2s kicks [see Eqs. 



( 24 ) and ( 26 ) with 6l) = 0] . We refer to this temporal recon- 
struction with period 2s di^ fractional antiresonant behaviour. 



In contrast, for even values of ^(1 - i), Eq. (33) reduces to 



Eq. (32), and the magnitude of v increases linearly with n at 
a rate proportional to 1/ ^[s. This leads io fractional resonant 
behaviour, where energy is transferred from the laser field to 
the atom cloud in quasi-periodic bursts. Fractional resonances 
are the main focus of this paper and are discussed in detail in 
Secs.|III]and[IVl 

As in our previous work (49) , we study the evolution of mo- 
mentum moments as a useful way to characterize the atomic 
dynamics. If the state of the atoms is initially prepared in the 
momentum eigenstate \k + P), the ^th-order momentum mo- 
ment after time nT is given by 

oo 

{p^^ = {nKf ^ \ckjiJ3, nT)Hj +I3f, (34) 
where co is determined by substituting y„/(Q,/5), for example 



as given by Eqs. ( 30 ), ( 3 1 ), ( 32 ), or ( 33 ), into Eq. ( 27 ). 



5. Evolution of an incoherent mixture 

To generalize the above treatment to the case of a cold 
thermal atomic cloud, we consider an incoherent mixture of 
plane waves with a Gaussian initial momentum distribution 
D{p) = DkQ3)/nK, i.e.. 



Dm 



1 



w y27T 



exp 



2w2 



(35) 



This corresponds to a Maxwell-Boltzmann distribution for 
free particles with temperature 7"^ = fP'K^w^ jMkB- 

The initial density operator for the system is taken to be 

/>-l/2 ^ 
dp\p)D{p){p\ = dpY \k+l5)Dkmk+P\. 
-'-1/2 ktZo 

(36) 

Evolving the density operator using Eq. (24), and considering 
the diagonal elements, we find that after time nT the momen- 
tum distribution evolves according to 1491 



AOS, t = nT)=Yj nTfDj{p\ (37) 



and the ^th-order momentum moment is given by 

>l/2 

i 

1/2 j,k=-o. 



Xl/2 «> 
dp y \c,j(Ji,nT)\ Dkmi+Pf 
1/2 



(38) 



6. Numerical implementation 

The numerical results presented in this paper are generated 
using a previously described Monte Carlo approach f49l. The 
initial condition consists of N plane waves distributed in mo- 
mentum according to D^ifi) of Eq. ( [35] ). The time evolution 
of each momentum eigenstate is evaluated by applying the 



Floquet operator ( 14 ), which is represented in a momentum 
basis 0. The final momentum distribution [see Eq. ([37])] is 
constructed by averaging the distributions resulting from each 
of the individual plane-wave time evolutions and collecting 
the data into a finite number of equal width momentum- space 
bins. The resolution of the averaged atomic momentum distri- 
bution can be improved by increasing the number of bins per 
unit momentum fiK. 



III. FRACTIONAL RESONANCES 
A. Fractional resonances in the p - subspace 

1. Fractional resonances with rational Q 

We first consider the dynamics of the atom-optical ^-kicked 
accelerator when the initial state of the system is a zero mo- 
mentum eigenstate, i.e., the dynamics are confined to the p = 
subspace. Figure [T] shows the resulting momentum distribu- 
tions, and corresponding second- and fourth-order momentum 
moments, when the pulse periodicity is equal to the Talbot 
time Tt = 2). We have chosen 0^ = O.Stt as an illustrative 
value typical of recent experiments [4,5,6,7,8,9,10,15, 24l . 

As a basis for comparison with the fractional quantum reso- 
nances, we show the known kicked-rotor behaviour of the Q = 
case [see Figs.[TJa) and[TJb)], i.e., a quantum resonance is 
observed and the ballistic expansion of the atom cloud is char- 
acterized by quadratic growth in the second order momentum 
moment. This growth is described by {p^}n = ff'K^cp^n^ /2, 
and is indicated by a solid line in Fig. [TJb) (see [49] and ref- 
erences therein). 

The quantum resonant dynamics are highly sensitive to the 
value of D, and, as predicted in Sec. |II B 2\ quantum an- 
tiresonance is observed for Q = 1 [see Figs. [TJi) and[T]j)]. 
Quantum antiresonance is characterized here by reconstruc- 
tion of the initial state every second kick, and the second- and 
fourth-order momentum moments oscillate accordingly (i.e.. 



between the ^ = and n = I values of Eqs. ( [39] ) and ( [40] ), 
respectively [49|). 

When the eff'ective gravity is chosen such that = l/s, 
where ^(1 - £) is even (see Sec. |IIB 4| ), we observe fractional 
resonances [see Figs, [ijc), [TJ^)' ^^^dJ^)]- ^ fractional 
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resonance, the energy transfer to the system is less efficient 
compared with the pure resonant case (Q = 0). The atomic 
momentum distribution expands in bursts of period Ht = s, 
and the quasi-periodic nature of the energy transfer is directly 
observable in the evolution of the second-order momentum 
moment [see Figs. [TJd), [TJf), and[TJh)]. The quasi-periodic 
behavior arises because, in the frame where the atom cloud is 
on average stationary, the phase of the optical standing-wave 
pulses varies from kick to kick with period Ht = s(l -\- sr mod 
2), as discussed in Sec. |IIB 3) 

For comparison, Figs.[TJk) and[TJl) correspond to an irra- 
tional value of Q. In this case, quasi-periodic behaviour is not 
observed. 



2. Momentum moment evolution 

For a zero momentum eigenstate, the dynamical evolution 
of the system is well characterized by the evolution of the 
second-order momentum moment, which, for Q = lis, even 
s{\ - £), and n an integer multiple of s, is given by |[63ll 



my 



2 ^y 

2s ' 



(39) 



Equation ( [39] ) describes quadratic growth in the second-order 
momentum moment (and therefore the kinetic energy) with a 
growth rate inversely proportional to ^.^ 

In the same parameter regime, the fourth-order momentum 
moment evolution for an initial zero momentum eigenstate is 
I63i 



30X 

8^2 



2s 



(40) 



To leading order, Eq. ( [4Q| ) describes quartic growth with a 
rate inversely proportional to s^. Therefore, the second- and 
fourth-roots of the second- and fourth-order momentum mo- 
ments, respectively, each grow (to leading order) linearly in n 
at a rate proportional to 1 / V^, as shown in Figs. [l|b), [l|d), 
[TJf), and[TJh). We observe that, when n is not an integer 
multiple of s, the momentum moments oscillate periodically 



around the analytic predictions ( [39] ) and ([40). 

For the infinitely narrow initial momentum distributions 
considered here, the evolution of the fourth-order momentum 
moment does not provide significant additional information 
over and above that extracted from the second-order momen- 
tum moment evolution (see Fig.[T]). However, as we show in 
Sec. |IIIB| the situation is quite diff'erent for initially broad, 
thermal momentum distributions, so we have included a de- 
scription of the fourth-order momentum moment dynamics 
here for completeness. 



^ For odd values of ^(1 - t) fractional antiresonances occur and all the mo- 
mentum moments return to their zero initial value every 2s kicks [see Eq. 
|33)]. 
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FIG. 1 : (Colour online) (Left-hand panel) Momentum distributions 
with a resolution of %K, and (right-hand panel) momentum moments 
of order {o) q - 1 and (x) q - A, for a (5-kicked accelerator. The 
initial condition is a zero-momentum eigenstate and parameters are 
M ^ \,T ^ Tt {i ^ 2), (pd ^ 0.8;r, and (a), (b) = 0; (c), (d) 
Q = 1/8; (e), (f) Q = 1/4; (g), (h) Q = 1/2; (i), (j) Q = 1; and (k), (1) 
Q = (1 + V5)/2. In the right-hand panel the solid and dashed lines 
correspond to Eqs. ( [39] ) and ( [40| ), respectively. The vertical lines in 
(d), (f), and (h), indicate where n is an integer multiple of s (as taken 
from Q = 1/5). 



3. Dependence on Q of resonant features 

In Fig. |2] we study the evolution of the second-order mo- 
mentum moment as Q is varied continuously from to 2 for 
both T = Tt/2 (£ = 1) and T = (^ = 2). The initial 
condition is in every case a zero-momentum eigenstate. We 
observe a rich structure of resonant and antiresonant features. 
Resonant (antiresonant) features appear dark (pale) and be- 
come narrower and more prominent with increasing n. The 
values of Q used in Fig. [T] are marked on Fig. |2jb), and the 
regions immediately surrounding these values are shown in 
more detail in Figs.[2jc) - (h). 

In particular, we see from Figs. |2ja) and|2jb) that chang- 
ing £ from 1 to 2 transposes the Q dependence of the various 
resonant and antiresonant features by one, as can also be de- 
duced from the structure of the transformed Floquet operator 
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FIG. 2: (Colour online) {p^)n^ in units of fiK for a (5-kicked accelerator initially in a zero-momentum eigenstate. Parameters are TV = 1, 
(pd - 0.8;r, and (a) T - Tj 12 {i - 1) and (b) T - Tj {i - 2). The markers highlight particular resonance features for T - Tj {i - 2), which 
are shown in more detail in (c) Q = 0, (d) Q = 1/8, (e) = 1/4, (f) Q = 1/2, (g) = 1, and (h) Q = (1 + V5)/2. 



([18]) [see Sec. |irB2l . 



B. Fractional resonances for a finite temperature cloud 

1. Manifestation of fractional resonances at finite temperature 

In the case of a finite-temperature cloud, where the ini- 
tial momentum distribution is described by the Gaussian ( [35] ), 
many quasimomentum subspaces are initially populated. As 
previously discussed for II = |49|, this results in some sig- 
nificant diff'erences in the overall finite temperature dynam- 
ics compared with the dynamics of a system restricted to the 

= subspace. 

Figure [3] shows, as a function of kick number, the momen- 
tum distributions and second- and fourth-order momentum 
moments for the parameters used in Fig.[T] with the exception 
that in Fig. [3] the initial condition in each case is a Gaussian 
momentum distribution with w = 2.5 (corresponding to a tem- 
perature of approximately 5 yuK in the case of caesium | 24 |). 
There are distinct diff'erences in this case compared to the zero 
momentum eigenstate evolutions shown in Fig.[T] For integer 
values of Q |[49j |63l, considered for the sake of comparison 
with the fractional resonant cases, pure resonant and antires- 
onant behavior occurs for the zero momentum eigenstate [see 
Figs.[TJa) and[TJi)]. At finite temperature the momentum dis- 
tribution behaves identically irrespective of whether D is even 



or odd [see Figs.[3ja) and[3];i)]. We observe that a small frac- 
tion of the cloud expands ballistically, but much of the distri- 
bution remains clustered near zero momentum. This behavior 
is due to resonances and antiresonances occuring concurrently 
in a finite temperature cloud, and has been well described pre- 
viously for D = [49]. 

For Q = 1/^, we again observe a small fraction of the cloud 
expanding ballistically at finite temperature [see Figs. [3jc), 
[3je),[3jg), and[3];i)]. However, as s increases, the rate of bal- 
listic expansion in the wings of the distribution is reduced. 
This is characteristic of fractional resonances and was more 
clearly observable in the zero-mometum eigenstate evolutions 
shown in Fig. [T] Also, with increasing s, the fraction of the 
cloud clustered near zero momentum becomes increasingly 
delocalized, as can be seen in Fig. [3] after 20 kicks. This in- 
dicates that, for rational values of Q, the antiresonances are 
not as eff'ective at localizing the atom cloud. A description of 
the fractional resonances at finite temperature is developed in 
SecIV] 

The evolutions of the second-order momentum moment are 
extremely similar for each value of illustrated in Fig. [3] We 
have found this to be a general property at finite temperature, 
which reflects the fact that resonant and antiresonant dynam- 
ics are occurring concurrently in the cloud. This rules out the 
second-order momentum moment as a useful means of con- 
cisely distinguishing difl'erent fractional resonances from one 
another, or indeed from the non-resonant case shown in Figs. 
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FIG. 3 : (Colour online) (Left-hand panel) Momentum distributions 
with a resolution of fiK, and (right-hand panel) momentum moments 
of order {o) q - 2 and (x) q - A, for a (5-kicked accelerator. The 
initial atomic momenta are distributed according to Eq. ([35]) and pa- 
rameters are N = 10000, w = 2.5, T ^ Tt {i ^ 2), 0j = 0.8;r, and 
(a), (b) Q = 0; (c),(d) Q = 1/8; (e), (f) Q = 1/4; (g), (h) Q = 1/2; 
(i), G) Q = 1; and (k), (1) Q = (1 + V5)/2. In the right-hand panel, 
the solid lines correspond to Eq. (|4TJ. The dashed lines in (b) and 
(j) correspond to Eq. ( [42] ), and the dashed line in (h) corresponds to 
Eq. ( [43] ). The vertical lines in (d), (f), and (h), indicate where n is an 
integer multiple of s (as taken from - lis). 



[3k) and|3];i). 

In contrast, the fourth-order momentum moment evolutions 
remain comparatively distinct, as shown in the right-hand 
panel of Fig. [3j We find that for larger values of s (where 
Q = 1/^), the fourth-order momentum moment increases at a 
slower rate. This is a quantitative signature of the fact that as 
s increases both the resonant and antiresonant dynamics be- 
come less efficient, and consequently the momentum distribu- 
tions become less peaked. The fourth-order momentum mo- 
ment evolution is sensitive to this and provides a useful way 
to characterize the manifestation of fractional resonances at fi- 
nite temperature. However, we note that the fourth-order mo- 
mentum moment evolutions for Q = [Fig.|3jb)] and Q = 1 
[Fig.[3]j)] are, like the second-order momentum moment evo- 
lutions, essentially indistinguishable. 



2. Evolution of momentum moments 

Deriving analytic expressions for the momentum moment 
evolution of a finite temperature gas is more involved than ar- 



riving at the zero temperature results (39) and (40). However, 
it is possible to derive such expressions for integer |63 1 and 
half-integer values of Q. In particular, we consider a thermal 
limit where w is taken to be large for the initial Gaussian mo- 
mentum distribution ( [35] ). 

For integer values of D, in the limit w » 1/ ^flni 1 631 , 



and 



(41) 



^ + + 3w^0> + 3w^ . (42) 



For half-integer values of (see Appendix [D]), in the limit 
w » 1 / 2 V27r^, and with the restriction that n is even, the evo- 
lution of the second-order momentum moment is again given 
by Eq. ( [4T] ), but the evolution of the fourth-order momentum 
moment is given by 
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(43) 

Consequently, and as observed in Fig.[3j the system dynamics 
for integer and half-integer Q values cannot be distinguished 
using {p^}n- However, the leading order term of the fourth- 
order momentum moment evolution is cubic in n, and the 
leading-order coefficient for half-integer Q is exactly half the 
leading-order coeflftcient in the integer D case. Therefore, it 
is possible to distinguish between fractional resonances with 
integer and half-integer values of D using {p^}n- 

Our numerical calculations indicate that, for rational D = 
r/s and in the large- w limit, {p^}n grows (on average) accord- 
ing to Eq. (41 ), and {p^}n grows to leading order (on average) 



cubically with n and with a growth rate proportional to 1/^. 
We emphasize that while we have only shown these trends 
to be exactly true for ^ = 1 and s = 2, the generalization 
to higher values of s is strongly supported by our numerical 
calculations. 



IV. EFFECT OF THE INITIAL ATOM CLOUD 
TEMPERATURE ON FRACTIONAL RESONANCES 

A. Momentum dependence of fractional resonance features 

1. Motivation 



In Sec. Ill A we considered the dynamics of the ^-kicked 
accelerator restricted to the JS = subspace, and observed 
fractional resonances for Q = 1/^ and T = Tt = 2). We 
found that the fractional resonances could be characterized 
by quadratic growth of the second-order momentum moment 
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at a rate inversely proportional to s, i.e., <^^)„ oc ri^/s [see 
Eq. ( [39] )] . Analogously, the fourth-order momentum moment 
evolves, to leading order, as oc n'^/s^ [see Eq. (40)]. 

In contrast to the zero-temperature limit, at finite tempera- 
ture (see Sec. |IIIB| ) we found that the second-order momen- 
tum moment grows linearly with ^, at a rate that appears to 
be essentially independent of Q, i.e, {p^}n oc n, and that the 
fourth-order momentum moment appears to evolve, to lead- 
ing order, cubically as oc ^^/^. 

The diff'erences observed between the zero-temperature 
limit and the thermal (large-w) limit arise because diff'erent 
quasimomentum subspaces evolve according to slightly dif- 
ferent (i.e., yS-dependent) transformed Floquet operators [see 
Eq. ([18])]. For a broad initial momentum distribution, all 



quasimomentum subspaces are populated and the observed 
momentum moment dynamics result from the appropriate av- 
erage over all the diff'erent /5-subspace evolutions [49] . It is 
therefore instructive to consider the dynamics of the diff'erent 
13 subspaces independently. 



2. Role of the quasimomentum (3 
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To better understand the yS-dependence of the system evo- 
lution, we investigate the evolution of momentum eigenstates 
\k + JS), where ^ = and e [-1/2, 1/2). In particular, we 
consider the evolution of the momentum variance 



FIG. 4: (Colour online) {{p^))]!^ in units of fiK for a (5-kicked ac- 
celerator where the system is initially prepared in the momentum 
eigenstate |/?). Parameters are A/" = 1, 0j = 0.8;r, T - Tj {i - 2), and 
(a) Q = 0, (b) Q = 1/8, (c) Q = 1/4, (d) Q = 1/2, (e) Q = 1, and (f) 

a = (1 + V5)/2. 



{{P^))n = {p\ - {P)l 



(44) 



For a momentum eigenstate, the variance displays clear sig- 
natures of the fractional resonances through the second-order 
momentum moment [see Sec. Ill A and Fig.[T|. However, the 
variance is more useful in this context because it evolves in- 
dependently of k [49], and its initial value is always zero irre- 
spective k and p. (The variance is introduced in a more general 
discussion of cumulants in Sec. |IVCT] ) 

Figure[4]shows the quasimomentum dependence of the evo- 
lution of {{p^))n^, for T = Tt = 2), and for a variety of dif- 
ferent values of Q. In Fig.|4ja), we show, for the sake of com- 
parison, the known results for the ^-kicked particle (D = 0): 
resonances occur foryS = and JS = ±1/2, and antiresonances 
occur for J3 = ±1/4 ||49l. For Q = 1 [see Fig.ge)], the reso- 
nance and antiresonance features are shifted in momentum by 
1 /4 compared to the Q = case. 

For = l/s, Figs. |4jb)-(e) show that increasing s in- 
creases the density in /3 of resonance and antiresonance fea- 
tures. However, the resonances and antiresonances remain 



evenly spaced and, as we discuss in Sec. IV A 4 the width 
of the fractional resonances are independent of s. 

In Fig.|4jf) we show the quasimomentum dependence of the 
evolution of {{p^))n^ for an irrational value of Q. We find that 
the resonance and antiresonance features observed for rational 
values become irregular and less well defined. 



3. Fractional resonance condition 

It is possible to quantify the p separation between fractional 
resonances. In the = subspace, for Q = l/^-, fractional 
resonances occur for even ^(1 - ^) (see Sees. |II B 4] and |III A| ). 
In general, for D = r/ it can be shown that fractional res- 
onances occur for even s[r - (1 + 2yS)^] 1631 . Inverting this, 
we find that fractional resonances occur for quasimomentum 
values 



~2r 



m 



(45) 



where m is integer. The fractional resonances are separated in 
momentum by 1 /is, breaking up momentum space into qual- 
itatively similar (although not identical) regions. The frac- 
tional resonances observed in Fig. |4] are consistent with Eq. 

<E3- 



4. Fractional resonance width 

Fractional resonances have a momentum width that de- 
pends on £ and the kick number n. As shown previously 
for the ^-kicked particle (491, an expression for the resonance 
width can be derived by considering antiresonance features in 
close proximity to the fractional resonances. Figure [5] shows 
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FIG. 5: (Colour online) {{p^))n^ in units of fiK for a (5-kicked ac- 
celerator where the system is initially prepared in the momentum 
eigenstate Parameters are A/" = 1, 0j = 0.8;7r, T - Tt {i - 2), 
and (a) Q = 0, (b) Q = 1/8, (c) Q = 1/4, and (d) Q = 1/2. 



a zoom in on theyS = fractional resonances from Figs.Qa)- 
(d). We observe hyperbolic curves of zero momentum vari- 
ance which we call reconstruction loci B9l . These are due 
to higher-order antiresonances which periodically reconstruct 
the initial state. 

The locations of the higher-order antiresonances that form 
the reconstruction loci can be determined for ^ = 1 and s = 



2 using Eqs. (30) and (31), respectively. The derivation is 
identical to that of the Q = case |49|, and we find that the 
reconstruction loci are described hy JS = mini, where m is an 
integer. The reconstruction loci closest to the resonance are 
then described by = ±1 Ini. 

We define the fractional resonance width as follows. The 
hyperbolic curves lying halfway between the fractional reso- 
nance and the reconstruction loci are the sequence of points 
at which the state begins the process of reconstruction. These 
transition loci are described hy p = ±l/2nL The resonance 
width can then be defined as the JS separation between the two 
transition loci adjacent to the fractional resonance, i.e.. 



FR 



1 

ni' 



(46) 



The fractional resonance width ( [46] ) is independent of s for 
s = I and s = 2. Our numerical calculations indicate that the 
width of the fractional resonances are, in general, independent 
of s, as illustrated in Fig. [5] 



B. Temperature dependence of the momentum moment 
evolution 

We investigate the eff'ect of the initial atom cloud temper- 
ature by considering the evolution of the second- and fourth- 
order momentum moments for Gaussian initial atomic mo- 



Figure |6] shows {p^}n and {p^}n for diff'erent values of w and 
Q. 

As w increases, for a given Q, the momentum moment evo- 
lutions make a gradual transition from the zero-temperature 



limit described in Sec. IIIB2 



limit behavior discussed in Sec. Ill A 2 towards the large-w 

The dashed lines in Fig. [6] 



indicate the zero-temperature limit, given by Eqs. (39) and 
( [40| ) for the second- and fourth-order momentum moments, 
respectively. The solid lines in the figure are lower bounds for 
the momentum moment evolutions in the large-w limit, i.e.. 



they are given by the large-w limit expressions in Sec. IIIB 2 
evaluated with w = 0. For the second-order momentum mo- 
ment, the solid lines correspond to Eq. ([4T] ) with w = 0. Note 
that for finite w, the large-w limit of (^^)^will be larger than 
that for w = 0, but will increase at the same rate. For the 
fourth-order momentum moment evolution, the lower bounds 
depend on Q and are defined by Eqs. ( [42| ) and ( [43] ), each with 
w = 0, for integer and half-integer Q, respectively. We have 
not evaluated analytic expressions for the lower bounds for 
general rational Q, although Figs.[6jf) and[6jh) strongly indi- 
cate that similar large-w limiting behavior exists. In addition 
to the general growth trends of the momentum moments, the 
quasi-periodic behavior with period = ^ is still observable 
(see Sec. [nB3] ). 

The deviation of the momentum moment evolution from the 
zero-temperature limit occurs at a particular kick number ^fr 
(indicated by the vertical dotted lines in Fig. [6]). As is well 
understood for the D = case (491, the energy transferred to 
the system becomes limited when the initial momentum width 
of the atom cloud becomes comparable with the momentum 
width of the resonance. If the initial momentum width is large 
compared to the resonance width, antiresonances play a role 
in the system dynamics and the momentum moment evolution 
tends towards the thermal large-w limit. In the case of the 
fractional resonances the same arguments apply, and indeed 
the width of the fractional resonances are independent of Q, as 
described in Sec. |IV A4l We define ^fr such that two standard 
deviations of the initial Gaussian momentum distribution lie 



within the resonance width ^y^FR (see Sec. IV A 4), i.e. 



^FR 



1 

4w^' 



(47) 



mentum distributions with diff'erent widths w [see Eq. ( 35 )] 



In Fig. [6] we also find that as s increases (where = l/s), 
the fractional resonances become increasingly sensitive to the 
initial temperature of the system, i.e., the momentum mo- 
ments approach their large-w limiting behavior for lower val- 
ues of w. This is most clearly observed in the second-order 
momentum moment, where the moment evolutions tend in- 
creasingly towards the solid line as s increases. We can 
explain this in terms of the number of resonance features 
spanned by the initial atomic momentum distribution. The 
fractional resonances are separated in momentum by 1/^^ (see 
Sec. |IV A"3] ). Thus, we can define a thermalization width Wth, 
above which the atomic momentum distribution initially spans 
more than a single resonance feature and the atomic dynamics 
reach a thermal limit. For = 1/s, this would be manifest 
in, for example, linear rather than quadratic growth in {p^}n 
with n, and cubic rather than quartic growth in {p^}n (see Sec. 
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FIG. 6: (Colour online) (Left-hand panel) {p^)n and (right-hand 
panel) for a Gaussian initial atomic momentum distribution 

with (X) w = 1/128, (o) w = 1/32, and (+) w = 1/8. Parameters 
are N = 10000, T ^ Tj {I ^ 2), 0j = 0.8;r, and (a), (b) Q = 0; 
(c), (d) a = 1/2; (e), (f) Q = 1/4; (g), (h) Q = 1/8; and (i), G) 
^1 = (1 + V5)/2. The dashed lines correspond to the w = analytic 
predictions [see (left-hand panel) Eq. ( [39] ) and (right-hand panel) Eq. 
([40|]. The solid lines correspond to the large-w limit lower bound 
[see (left-hand panel) Eq. (|4lJ, (b) Eq. (|42}, and (d) Eq. ([43]), all with 
w = 0]. The vertical dotted lines indicate n - n^^ [see Eq. ([47|], 
and the solid vertical lines in (c)-(h) indicate where n is an integer 
multiple of s (as taken from Q = 1 /5). 



IIIB 2| ). Quantitatively, we define Wth such that one standard 
deviation of the initial Gaussian momentum distribution lies 
within the fractional resonance separation, i.e.. 



C. Momentum cumulant evolution 



7. Motivation and background 



In Sec. IIVBI we found that the evolution of the momen- 



(42), 



tum moments depends explicitly on w [see Eqs. (41 
and ([43])]. For this reason, we defined a lower bound for the 
large- w limit of the momentum moment evolutions. As can be 
shown explicitly for integer 1631 and half-integer values of Q, 
this difficulty does not arise if we consider the change in the 
momentum cumulants from their initial value. In this way it 
is possible to have a well defined high-temperature limit. 

A ^th-order momentum moment {p^}n is dependent upon 
all moments up to order (^ - 1). In particular, {p^}n and (p^)^ 
are not independent quantities. Using an iterative process, mu- 
tually independent cumulants {{p^))n can be constructed from 
the moments {p^)n, . . . , |58, 65, 66, 67 1. The first- 
order cumulant {{p))n is simply the mean {p)n. The second- 
order cumulant {{p^))n is the variance [as defined in Eq. (44 )]. 
The fourth-order momentum cumulant, the kurtosis, is 

«/» = {p')-A{p'){p) + l2{p^){pf -?>{p^f -6{p)\ (49) 



where we have dropped the n subscripts for brevity. 

For a symmetric momentum distribution, the odd-ordered 
momentum moments are identically zero, and the second- and 
fourth-order cumulants simplify to 



{{P'))n ={p')n, 

{{p'))n ={p\ - 3{p^}l 



(50) 
(51) 



If the momentum distribution is symmetric initially, one can 
show that, for = r/ s, it must evolve to a symmetric distribu- 
tion at intervals of s kicks [63 1. One can therefore infer (and 
indeed observes) that the odd moments should be bounded in 
value, and remain relatively insignificant in comparison to the 
rapidly growing even moments. We therefore, at times, make 
use of Eq. ( [50] ) and Eq. ( [5T] ) even when their use is not fully 
justified. Numerically, the fact that we employ a Monte-Carlo 
method yields small, although in general non-zero, odd mo- 
ments. Therefore, in numerical calculations we evaluate the 



momentum moments using Eq. (38), and then determine the 
cumulants from the moments using Eqs. ([44]) and ([49]). 



1 

2fs' 



(48) 



Equ ation ([48] ) is consistent with the large-w limit described in 
Sec. Ill B 2| (and taken more formally in Appendix [P]), where 



terms involving exp(-7r^w^/2w^j^) were taken to be negligible. 

Finally, we note that if Q is chosen to be irrational [see Figs. 
[6ji) and|6]j)], we observe that the second-order momentum 
moment grows linearly with n, and the fourth-order momen- 
tum moment grows quadratically. The data becomes smoother 
as w increases, but the growth rate appears to be largely inde- 
pendent of w. 



2. Zero -temperature limit 

Before considering the high-temperature limit of the cu- 
mulant dynamics, it is first instructive to consider the zero- 
temperature limit where the initial state is a zero-momentum 
eigenstate. The initial state is symmetric in p so {{p^))n = 
{p^)n is given by Eq. ( 39 ) which is valid for n an integer mul- 



tiple of s (where Q = 1/^). Substituting Eqs. ([39]) and ([40) 
into Eq. ( [5T] ) yields the known result |63 1 



8^2 



2s 



(52) 



where again n is an integer multiple of s. 
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FIG. 7: (Color online) {{p"^))]!^ for a (5-kicked accelerator with Gaus- 
sian initial momentum distribution of standard deviation {o)w - 1/4, 
(□) w = 1/16, (X) w = 1/32, ( ) w = 1/64, (a) w = 1/128, and (+) 
w = 1/1024. Parameters are N = 10000, T = Tt = 2), 0j = O.Stt, 
and (a) Q = 0, (b) Q = 1/2, (c) Q = 1/4, (d) = 1/8, and (e) 
Q = (I + V5)/2. The dashed lines correspond to Eq. \52\ , and the 
solid lines correspond to (a) Eq. ( [54| and (b) Eq. ([55]). The vertical 
lines in (b)-(d) indicate where n is an integer multiple of s (as taken 
from^l = 1/^). 



3. High-temperature limit 

For a Gaussian (and therefore symmetric) initial mo- 
mentum distribution of standard deviation w, by definition 
{{P^))q = (fiRfw^ [see Eqs. Jsol and (|E21]. In the thermal 



limit, where w is taken to be large, the momentum variance is 
given by Eq. ([4T]). Subtracting the initial value gives 



(53) 



which is independent of w. Note that Eq. ( |53) ) applies to both 
integer and half-integer values of Q (although for half-integer 
values of Q it applies only for even values of n). Our numer- 
ical calculations strongly indicate that Eq. ([53] ) should apply 
for all values of Q (see Fig.[3]and Sec. IIIB 2[ ). 

A Gaussian distribution has no non-zero cumulants of order 



greater than two so, when considering an initially Gaussian 
momentum distribution, we have {{p^}}o = 0. For integer 



values of D, we substitute Eqs. (41 ) and (42 ) into Eq. ( 51 ) and 
find that 



1= 



3^> 



n + 



(54) 



Equation ( |54] ) is also independent of w. An analogous expres- 
sion for half-integer values of Q can be found by substituting 
Eqs. ([41]) and ([43]) into Eq. ([51]) to yield 



1'- 



n + 



(55) 



which is valid for even n. Equations ([54]) and ( [55] ) show that 
the cubic leading-order cumulant behavior, and its coefficient, 
is unchanged compared to the corresponding momentum mo- 
ments [see Eqs. ( [42] ) and ([43])]. For general rational values 
of II, we expect that the change in the momentum cumulants, 
compared to their initial values, will again be independent of 
w, for w sufficiently large. 



4. Temperature dependence of the momentum cumulant evolution 

The second-order momentum cumulant evolution is essen- 
tially described by the second-order momentum moment evo- 



lution [see Eq. ( [44] ) with negligible {p)n\- Therefore, the tem- 
perature dependence of the second-order momentum cumu- 
lant evolution does not provide any further information than 
that presented for the second-order momentum moment evo- 
lution in Sec. lIVBl In this section we concentrate on the tem- 
perature dependence of the fourth-order momentum cumulant 
evolution. 

Figure [T] shows the evolution of the cube-root of the fourth- 
order momentum cumulant for diff'erent values of w and Q. 
The fourth-order momentum cumulant initially evolves ac- 
cording to the zero-temperature limit, indicated by the dashed 
lines in the figure, and is negative indicating a broad momen- 
tum distribution characteristic of quantum resonance phenom- 



ena. At the kick number n = ^fr [see Eq. (47 )], the cumulants 
deviate from the zero-temperature limit and eventually be- 
come positive, indicating a sharply peaked momentum distri- 
bution which is caused by antiresonances playing a significant 
role in the dynamics. For irrational values of D [Figs. [7];e)], 
the absence of resonance and antiresonance features means 
that the fourth-order cumulant remains negative. 

For sufficiently large values of w, the fourth-order momen- 
tum cumulant tends to the thermal limit, indicated by the solid 
lines in Figs.[7ja) and[7];b). In the thermal limit, the fourth- 
order momentum cumulant evolution is independent of w and 
is characterized by cubic growth with n [see Eqs. ([54]) and 



( 55 )]. The growth rate is smaller for higher values of s (where 
= 1/s), indicating the less peaked momentum distributions 
of higher-ordered fractional resonances (see Fig. [3]). 

We observe that for larger values of s (where = l/s), 
the cumulant evolution tends to the large-w thermal limit be- 
havior more readily, i.e., for smaller values of w. We have 
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FIG. 8: (Color online) The asymptote gradient bopt of ((p^)y^^. Pa- 
rameters are N = 10000, T = Tt (i = 2), 0j = O.Stt, and (x) Q = 0, 
( ) Q = 1/2, (o) Q = 1/4, (□) Q = 1/8, and (a) Q = 1/16. The 
dashed lines are shown to guide the eye. The horizontal solid lines 
correspond to the thermal asymptote gradients of Eqs. ( [54| and ([55]). 
The vertical dotted lines correspond to Eq. ([48]). 



investigated this quantitatively by determining the asymptotic 
dependence of the fourth-order cumulant. This was done by 
fitting a particular number of asymptotic points using a least- 
squares method according to {{p^}}n^ = bn-\-c. The most ap- 
propriate asymptote, with b = /?opt, was identified by choosing 
the number of asymptotic points such that the standard error 
in the gradient b was minimized. 

In Fig. [8] we show the optimum asymptote gradient /?opt as 
the width w of the Gaussian distribution is varied. We find 
that as w increases, the asymptote gradient saturates to its 
thermal limit value /?th- For integer and half-integer values 
of Q, bth can be determined from Eqs. ( [54) and ( [55] ) to be 
hh = {^yAY'\nK)^l^ and /?th = (04/8)^(;i^)4/^7espec- 
tively, as indicated by the solid horizontal lines in Fig. [8] 

As previously indicated, the asymptote gradient saturates 
to the thermal value /?th more readily for higher-order frac- 
tional resonances, i.e., for larger values of s (where = l/s). 
This saturation occurs when the atomic momentum distribu- 
tion initially spans a single fractional resonance feature in j3, 
as described in Sec. |IVB[ Quantitatively, the thermalization 
width Wth is defined by Eq. ( [48] ), and as indicated by the verti- 
cal dotted lines in Fig. [8] it accurately predicts the temperature 
at which the fractional resonant dynamics saturates in the ther- 
mal limit. 



V. DISCUSSION AND CONCLUSIONS 



eff'ect of the initial atom cloud temperature on fractional quan- 
tum resonant dynamics, and have shown that the kick num- 
ber at which the system dynamics clearly deviate from the 
zero-temperature case is identical for all fractional resonances 
and is inversely proportional to the initial atomic momentum 
width of the cloud. However, higher-order fractional reso- 
nances are more sensitive to thermal eff'ects and the dynamics 
saturates to the thermal limit behavior at lower temperatures. 

The observed resonant features are sensitively dependent 
on the diff'erence between the local gravitational acceleration 
and the applied walking- wave acceleration. As calibration 
of the phase shift between the two counterpropagating waves 
making up the walking wave to between parts-per-million and 
parts-per-billion is in principle possible 1681 . there may be a 
potential application in the sensitive atom-optical determina- 
tion of the local gravitational acceleration |9, 5X 68^..69J. For 
a pulsing periodicity of equal to the Talbot time, one would 
for example observe resonant behaviour when the local grav- 
itational acceleration is exactly cancelled by the acceleration 
of the walking wave. 

Typically in cold atom experiments the initial atomic mo- 
mentum width is on the order of w ~ 0.04 for Bose-condensed 
rubidium |25|, and w ~ 0.008 for Bose-condensed cesium 
155 J . At these temperatures the experimental observation of 
fractional resonances is accessible. 
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APPENDIX A: FACTORIZATION OF THE 
FREE-EVOLUTION OPERATOR 

To factorize the free-evolution component of the Floquet 
operator ([9]) into position- and momentum-dependent parts, 
we first consider a general operator of the form 

d(T) = exp [prp^ -h cTT^p -h T]T^^ Qxpiyrz). (Al) 

Taking the partial derivative of O with respect to r gives 



dO(T) 



-h [pp^ -h ICTTp -h St/T^)] 0(t). 

Using the general operator identity 1701 



We have presented a detailed theoretical investigation into 
fractional quantum resonances in the atom-optical ^-kicked 
accelerator. Fractional resonances occur for particular ratio- 
nal values of the scaled eff'ective gravitational acceleration Q, 
and are characterized at zero temperature by expansion of the 
atom cloud in quasi-periodic bursts. We have considered the 



e^Be-^ =B + [A,B] + — [A, [A, ^]] -h • • • (A3) 

(truncated to the first two terms, as all further terms are zero), 
we find that 

(A4) 
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Hence, ( |A2| ) simplifies to 
dd(T) 



dr 



= [^p^^ + 2t (cr - ifipy) p + yz + T^ Or] - ificry)^ d(T). 

(A5) 

Choosing cr and rj such that terms proportional to r and 
vanish, i.e., cr - ifipy, and subsequently 7/ = -ff-py^/S, Eq. 
( |A5| ) reduces to 



— ^ - (p^2 + +rz)0(T). 

or 



(A6) 



Noting from Eq. ( |A1| ) that 0(0) = 1, Eq. (|A6|) can be inte- 
grated to give 

6(t) = Qxpiprp^ + yrz). (A7) 

Hence, 

exp (^prp^ + yrz) = exp [prp^ + ifipyr^p - ff'py^r^ 



X exp (yrz) . 



(A8) 



Setting p = -ijlMfi, y = -iMa/H, and r = T, yields the fac- 
torized form of the free-evolution component of the Floquet 
operator, as used in ([TT]). 



APPENDIX B: ^stat FOR RATIONAL a 

For particular values of Q, Hamiltonian ([19]) has a period- 
icity in addition to the delta-kicks. This arises due to the Q- 
dependent phase of the standing- wave pulses, which is identi- 
cal for a given remainder of n^Q./2. In the case where = r/s, 
this occurs whenever rPr/s mod 2 takes the same value. 

First we consider the case where sr is even. At the ^th kick, 
i.e., where n = s,wq have that ri^r/s mod 2 = sr mod 2 = 0, 
and the phase of the standing-wave pulse is identical to the 
phase taken for ^ = 0. For n = s -\- n\ 

(s-^n'f- mod 2 = (2sn' -^n'^)- mod 2 = n'^- mod 2, (Bl) 



Theefore, we must calculate v of Eq. ( [28] ). In particular cases, 
the Gauss sum CSJ [73l |74j| in y can be evaluated analyt- 
ically t63\ . Here we evaluate v for half-integer values of Q, 
i.e., Q = r2/2. In this case 



n-l 



V2/2^^20T^-/;r/r2/2^ (CI) 



7-0 



For even values of j, exp(-/7r/r2/2) = 1, and for odd values 
of J, exp(-i7rfr2/2) = exp(-/7rr2/2). For an even number of 



kicks n = 2n', Eq. (CI ) can be rewritten as two sums, each 
with the same number of elements, i.e.. 



Ain'y 2i7m'^r2 



n'-l 
71=0 



,4/71 T 



n'-l 

-inr2/2 ^ ^2/(272 + l)T 
72=0 



(C2) 



Simplifying, we find that 



yn-2n'Ar2/2,/3) = /."^^"^'^ [l + .^(^^'-Z^)] e'^^\ (C3) 



7=0 



Evaluating the geometric sum, and simplifying further, yields 

sin(2^'T) (C4) 
xcos(T-7rr2/4) — 



sin(2T) 



Resubstituting for n = 2n' gives 



ynAr2/2,fi) = 2/.-"t^^"^i>^-^/4] _ ^^^^^^ ^A^'^nT) 



sin(2T)' 
(C5) 



which applies for even n only. 



and, in the frame where the cloud of atoms is (on average) sta- 
tionary, there is an exact ^-periodicity to the time dependence 
of the spatial phase of the standing- wave pulses. 

Similarly, if sr is odd, at the 2^th kick, ri^r/s mod 2 = 
4sr mod 2 = 0. For n = 2s -\- n\ 



(2s -h n'f- mod 2 = n'^- mod 2, 



(B2) 



and there is a 2 ^-periodicity to the spatial phase of the 
standing- wave pulses. 

A unified expression for the temperal period ht of the phase 
of the standing- wave pulses (valid for even and odd sr), is 
given by = '^'(1 + sr mod 2). 



APPENDIX C: EVALUATING v FOR HALF-INTEGER Q 



To determine the time evolution of the momentum eigen- 



APPENDIX D: MOMENTUM MOMENTS FOR 
HALF-INTEGER Q 

At finite temperature, the ^th-order momentum moment is 



given by Eq. (38). Inserting the matrix elements from Eq. 



state ( 21 ) requires knowledge of the complex amplitudes ([26] 



( [26] ), using that oj = (pd\y\ [see Eq. ([27])], and returning to the 
atomic momentum notation p = fiK(k -h yS), we find that 

/oo 
dp J] J] i<l>d\y\) Dip) mj + pf . (Di) 



Here we consider only the second- and fourth-order mo- 
mentum moments in detail. The second-order momentum mo- 
ment is given by 

/ oo 
dp Yj J] ^^d\y\) Dip) [inKff + IHKjp + /] . 
j= — oo 

(D2) 
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Using the Bessel function identities Ti'JL-oo j J ji^) = 7/2 



I.7=-^jj]{r]) = 0, and Y,%_^J]{r]) = 1 |49l|64l, Eq. (|D2) 
becomes 

{p')n = j dpD{p)p^ + (nKf'^ j dpD(p)\v\\ (D3) 
The fourth-order momentum moment is given by 

dp y jj(4>Ay\)D(p)\{nK)'f 

j~i (D4) 

+A(nKffp + GihKffp^ + 4hKjp^ + /] . 

Using the Bessel function identities above, as well as 
I.7=-^fj]iri) = and i:%_^fjj(il) = 3;/V8 + 7]^/2 m, 
Eq. (|D4|l becomes 



{p\= j p'D(p)dp + (nK)''^ I dpD(pM 



(D5) 



+ {nKf'icpl j dpD{p)pM^ 

+ {nKf^-^ j dpD(p)\v\\ 

For half-integer values of Q, the second-order momentum 
moment evolution for a Gaussian distribution is given by Eq. 



( D3| ) with V of Eq. ([31). Using Eqs. (|E2]l and ( |ET0| >, we find 
that 



j;2 n/2-1 



(D6) 



m=l 



In the Hmit w » 1 / 2 V27r^, the decaying terms become negU- 
gible and we can write 



2 



Note that Eqs. ( |D6| ) and ( |D7| ) only apply for even n. 

The fourth-order momentum moment for a Gaussian distri- 
buti on is given by Eq. (|D5]). Using Eqs. (|E3]), ( |E10| ), ( |E14| ), 
and ( |ET9l ), we find that 



=3h.4 + ^„ + 3vv2^2„ + ^„(„2 + 2) (3 - 

«/2-l 



.2 ...2/72,,, 2x^-87:2^2^2^2 



1 



=1 



- nm^ - m-\ — n(n^ -\- 2) 
6 



X ^^^-^^^f^^^^^^ - ^-^7T^(m+lfi^w^ _ ^-8;r2(m-l)2^2^2j 
Q j4 n-2 r ^ 

30^ ^ w 



^ -h ^(3^ - l)m -(1 - ^ ) 



|^6^-8;r2m2 



^2>v2 _ ^-8;r2(m+l)^rw2 _ ^-8;r2(m- 



- 1)2^2^2 j 



(D8) 



In the limit w » 1/2 V27r^, we can write 

{P'^)n ^ 4 , , o ,2 2 , / 2 , c\ rT^n\ 



Again, we emphasize that Eqs. ( |D7| ) and ( |D9| ) only apply for 
even n. 



APPENDIX E: INTEGRALS OVER D{p) 

1. Initial momentum moments 

At finite temperature, the initial ^th-order momentum mo- 
ment is given by Eq. ( |D1| ) with y = 0, i.e., 

{p^}o = J dpD(p)p^. (El) 

In the case where the atomic momenta are distributed accord- 



ing to the Gaussian (35 ), the initial second-order momentum 
moment is 

{p\ = (nKfw\ (E2) 
and the initial fourth-order momentum moment is 

{p\ = XnKfw\ (E3) 

2. Integrals involving y 

a. Evaluating J dpD{p)\v^ 
When takes half-integer values, i.e., Q = r2/2 with r2 



odd, y is given by Eq. (31 ) for even kick numbers n = In' . 
Hence, 



r o r sin^(2^'T) . 

(D7) dpD{p)\v\^=A dpD(p) J -cos\r-7ir2/4). 

J J sin (2T) 



Using the expansion f63| 
sin^(^x) 



(E4) 



= n + 2^^(n - m) cos(2mx), (E5) 



sin (x) 

cos(4mT) = cos(4m7i£p/fiK), and 

cos^(Y-^Ui + i(-l/^<->'/^ sinful 



4 1 2^2^-'^ 



\ nK 



(E6) 



equation (E4) becomes 



r dpDipW^ =2n'+4"^(n'-m) f dpD(p) cos I ^^^] 
+ 2(-l/-('--i)/2j apD(p)sin{^^^ 

X -h2^(fz' -m)cos|^^^^^| , 



(E7) 
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where we have used the normalization of D(p) to evaluate the 
first term. In the case where the atomic momenta are dis- 



together with 



tributed according to the Gaussian (35), the first integral on 



the right-hand side of Eq. ( E7 ) can be evaluated using 

dpcos{Ap)e-P''^' = ^J^ae-^'^''\ (E8) 



cos 



3 1 



8 ^\~nK] 



■ cos 



The remaining integrals vanish due to the even parity of D(p). 
Therefore, we find that 



^+(r2-l)/2 



sm 



27r£p 
fiK 



(E17) 



n'-l 



r dpD{p)\v\^ = 2^' + 4 - m)e-^^'^'^'''' 

^ m=l 

Substituting n = 2n' yields 

p n/2-l 

j dpD(p)\y\^ = ^ + 2 2 - 2m)^-^'^'^'^'^' 
where n is even. 



m=l 



(E9) 



(ElO) 



and following the treatment in App. E2a we find that 



dpD(p)\yf =-n\2n'^ + 1) [3 - ^"^^'^'^'j 

2n' 



b. Evaluating J dpD(p)p^\y\ 
For Q = r2/2 and n = 2n\ 
fdpD(p)p'\rf'4f 



- 2n'm^ - m + ^(2^ + 1) 



dpD(p)p^\yf = 4 I dpD(p)p^ ^^^ cos^(T-7rr2/4). 
' sin2(2T) 

(Ell) 

Following the treatment in App. |E2a and using (E2 ) and the 
integral 



e 

In'-l 



J / cos(Ap)e-P'^'''dp = i V^Qf^ |l 



we find that for Gaussian D(p) [see Eq. ([35])], 

r 

2n' -\-4^(n' -m) 



J dpD(p)p^\yf =(nK)W 



m-l 



(E12) 



(E13) 



X (1 - 4mV^V2)e-^'"''^^'^' 



Substituting n = 2n' yields 



dpD(p)p^\y\^ =(nK)W 



n/2-l 



n + 2 ^ (n - 2m) 



(E14) 



-2 



3 



— -2n'n? + -(12^'2 - l)m 



[6^ 



+1^(1 -4«'2) 



(E18) 



Substituting ^ = 2^' then yields 



n/2-1 



- ^m^ - m + -n{r^ + 2) 
6 



[6e- 



-^n'^im-lfeV 



c. Evaluating J dpD(p)\y\'^ 
For Q = r2/2 and ^ = 2n\ 



J dpD(p)\yf = 16 J dpD(p) 



sin^(2^^T) 
sin'^(2T) 



cos'^CT - 7rr2/4). 

(E15) 



Using the expansion |163 



sin (nx) n 2 



n-l 



sin^(x) 3 



(2n^ + l) + ^cos(2mx) 
m -2nm -m+—{2n + 1) 

2«-2 



^ cos(2mx) 



— - 2nm^ + 14^^ - - Im + y (1 - 4n^) 



(E16) 



«-2 
m=n/2 



.3(1-.^) 



_ ^-&K^(m-H)-(-w- _ ^-&K^im-l)-(-M^ 



(E19) 
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